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ABSTRACT 



I present a simple calculation of the expected mean CO brightness tempera- 
ture from the large scale distribution of galaxies during cosmic reionization. The 
calculation is based on the cosmic star formation rate density required to reion- 
ize, and keep ionized, the intergalactic medium, and uses standard relationships 
between star formation rate, IR luminosity, and CO luminosity derived for star 
forming galaxies over a wide range in redshift. I find that the mean CO bright- 
ness temperature resulting from the galaxies that could reionize the Universe at 
^ = 8 is Tb ~ l.l(C/5)(/esc/0.1)~^ /iK, where fesc is the escape fraction of ioniz- 
ing photons from the first galaxies, and C is the IGM clumping factor. Intensity 
mapping of the CO emission from the large scale structure of the star forming 
galaxies during cosmic reionization on scales of order 10^ to 10^ deg^, in combi- 
nation with HI 21cm imaging of the neutral IGM, will provide a comprehensive 
study of the earliest epoch of galaxy formation. 

Subject headings: cosmology: theory, reionization - galaxies: formation - galaxies 



1. Introduction 



Cosmic reionization corresponds to the epoch when the neutral intergalactic medium 
(IGM) is reionized by light from the first galaxies (Loeb & Barkana 2001). This epoch, and 
the preceding 'dark ages' between recombination and reionization, represent the 'Universum 
Incognito' - the last unexplored epoch of cosmic structure formation. Current observational 
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constraints, based primarily on large scale polarization of the CMB and Lya scattering of 
light from z > Q quasars (the Gunn- Peterson effect), suggest that rcionization may have 
had a large variance in space and time, starting as early as ^ ~ 15, and extending down to 
z ~ 7 (Fan, Carilli, Keating 2006). Interestingly, a recent census of star formation in the 
first galaxies suggests that, if rcionization is due to the galaxies that are seen in current deep 
near-lR fields (with modest extrapolation to lower luminosity; eg. Bouwens et al. 2010), the 
dominant epoch of rcionization could occur at z < 9 (Robertson et al. 2010). 

The most direct method for studying the evolution of the neutral IGM is through the 
21cm line of neutral hydrogen. Extensive theoretical and observational work is on-going 
to study the HI 21cm signal from the neutral IGM during cosmic rcionization, into the 
preceding dark ages (Furlanetto et al. 2006; Morales & Wyithe 2010). 

An equally important, and complementary, area of study will be to trace-out the large 
scale structure of the star forming galaxies that reionize the Universe. Indeed, a cross 
correlation between the large scale structure in star forming galaxies, and the structure of 
the neutral IGM, will provide a comprehensive view of Universal evolution during the epoch 
of first light. Relative to each other, these two probes provide 'inverse views' of the Universe 
(Lidz et al. 2009). 

The difficulty in determining the galaxy distribution is the very large scales involved. 
The HI 21cm studies are probing scales ranging from 100 deg^ to 1000 deg^, with depths of 
z ~ 6 to 10 (Parsons et al. 2009; Lonsdale et al. 2009, Harker et al. 2007). Covering such 
large volumes with near-IR surveys of 2; > 7 galaxies is well beyond any instrumentation in 
the foreseeable future. Moreover, studies of rest-frame UV light and/or the Lya emission 
from galaxies well into rcionization likely provide a biased view of the galaxy distribution 
due to the combination of strong Lyo; absorption by the neutral IGM, and any dust in the 
host galaxies (Fontana et al. 2010). 

In this letter, I consider the possibility of tracing the large scale structure of star forming 
galaxies during rcionization using the technique of intensity mapping of the molecular gas 
in star forming galaxies, as traced via the CO emission lines. Intensity mapping involves 
imaging the aggregate emission from thousands (or more) of galaxies on very large scales 
(hundreds of comoving Mpc). Individual galaxies are not detected, just the summed signal 
from the large scale distribution of galaxies. This technique has been explored using the HI 
21cm line out to 2; ~ 1 (Chang et al. 2010). Note that the ALMA and the EVLA are able to 
detect the CO emission from individual, massive galaxies beyond z ~ 6 (Wang et al. 2010; 
Carilli 2010; Walter & Carilli 2007), however, the field of view of these instruments is small 
(< 1'), precluding surveys on scales of tens of degrees. Intensity mapping offers a straight- 
forward alternative for tracing the very large structure in the galaxies that reionize the 
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neutral IGM, using smaller area telescopes, without the need to detect galaxies individually. 

My goal is to estimate the mean brightness of the CO emission on large scales from 
the galaxies that reionize the neutral IGM. The notion that molecular line emission from 
galaxies at z = to 10, could be a substantial CMB foreground was considered by Righi et 
al. (2008) using an extended Press-Schecter formalism for merging halos to set the evolution 

of galaxies, with analytic recipes for the implied CO emission from these galaxies. In this 
letter, I take a simpler approach, based on the required cosmic star formation rate density to 
reionize the Universe. This approach allows for a straight- forward derivation of the mean CO 
surface brightness from the galaxies that could reionize the neutral IGM at a given rcdshift, 
although it does not predict the structure in this emission. Clearly, the uncertainties in the 
assumptions are significant, and these are pointed out below. However, the estimate is meant 
to be order of magnitude, and useful for initial consideration of future instrumentation to 
perform CO intensity mapping. 

2. Calculation 

Following is a calculation of the mean CO brightness expected from the galaxies that 
could reionize the Universe at a given redshift. I adopt the standard concordance cosmology 
(Spergel et al. 2007), and comoving coordinates are used throughout. The basic information 
that goes into the calculation can be summarize as follows: 

• The cosmic star formation rate density required to reionize (and keep ionized) the 
neutral IGM. 

• A conversion from star formation rate to IR luminosity based on known properties of 
galaxies. 

• A conversion from IR luminosity to CO luminosity based on known properties of galax- 
ies, at least out to z ^ 3. 

2.1. Star formation rate density and IR luminosities 

Madau, Haardt, and Rees (1999) calculate the star formation rate density {psFR in 
Mq yr~^ Mpc~^, comoving) required to reionize the IGM. The basic assumptions are the 
mean baryon density, the clumping factor (recombinations), the escape fraction of ionizing 
radiation from galaxies, and a Salpeter IMF over 0.1 to 100 Mq range. The Madau et al. 
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calculation has been updated to the latest cosmological parameters by Bunker et al. (2010) 
in their section 3.3. This relation is key to the calculation below, and I repeat it here: 

(1) = SFR/Volume = 0.005 /"J [^^^]^ {C/b)X yr"^ Mpc'^ 

where /g^e is the ionizing photon escape fraction, C is the clumping factor, and X is a factor 
that depends on the cosmic baryon density and h^, and equals 1 for the standard cosmology 
(X will be removed from hereon). Simulations suggest a value of C somewhere between 5 
and 30 (Furlanetto et al. 2006). The value of f^sc ~ 0.06 for the Milky Way and other nearby 
star forming galaxies (Putman et al. 2003). Shaplcy et al. (2006) show that this may rise 
to 0.1 to 0.2 for 2; ~ 3 LBGs (see also Nestor et al. 2011), and there is marginal evidence 
that this may increase further to higher redshift (Bouwens et al. 2009). 

The relationship between far IR luminosity (FIR) and star formation rate has been 
considered in numerous studies over the last two decades. I adopt the standard relationship 
for nearby galaxies given in Kennicutt (1998): 

(2) Lfir = 1.1 X 10^° SFR Lq 

where SFR is in M© yr~^. I will assume in the following calculation that this relationship 
applies on average to all galaxies, and hence applies to the volumetric average over galaxies 
on large scales. 

In this case, the SFR as a function of LpiR can be solved for in equation (2), and 
substituted into equation (1), yielding the FIR luminosity per unit cosmic volume required 
to reionize the neutral IGM: 

(3) pFin = 5.5 X 10' f'l [^^]' iC/5) MpC' 

where p in this, and subsequent, equations indicates the comoving volume density of the 
subscripted quantity, in this case, the FIR luminosity. 

2.2. Cosmological relationships for CO and FIR luminosities 

The 'integrated Kennicutt- Schmidt law' relates the FIR luminosity of galaxies to the 
CO luminosity, and also has been a topic of extensive study for both nearby and distant 
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galaxies. The underlying physical relationship is between the star formation rate and the 
gas mass, but for the purposes of this letter, the critical component is the observational 
relationship. 

In this paper, I employ the latest empirical relationships between CO and FIR luminosity 
derived by Daddi et al. (2010) for both nearby and distant star forming galaxies. They 
conclude that there are two, roughly linear, relationships between CO luminosity and FIR 
luminosity for galaxies, one relevant for normal star forming galaxies at both low and high 
redshift, including the sBzK galaxies at z ^ 2, and one relevant for dense starburst galaxies, 
such as the ULIRGs nearby, and the submm galaxies and quasar hosts in the distant Universe. 
I will use the median ratio they derive for the dominant, gas rich, star forming disk galaxy 
population at low and high redshift: 

(4) L'co = 0.02Lfir K km s~^ pc^ 

where Lpm is in Lq. 

The units for L'qq were originally designed for spatially resolving observations of CO 
in the Galaxy, where brightness temperature was paramount, and from which the empirical 
relationships between molecular gas mass (H2) and CO luminosity have been calibrated. I 
return to this point below. 

We next consider the standard relationships between CO luminosity and observed flux 
from Solomon & Vanden Bout (2006). Since these relations are key to the calculation, I 
repeat the critical relations here. 

Solomon & Vanden Bout relate CO luminosity to flux density and line width as follows: 

(5) Leo = 1-0 X 1035Al^(l + z)-^VrDl Lq 

where the luminosity distance, Dl, is in Gpc, Vj. is the rest frequency in GHz, flux density, 
S, is in Jy, and velocity width, AF, is in km s~^. They also give a related equation for L'^q-. 

(6) L'co = 3-3 X 10^^SAVDliy;\l + z)'^ K km s"^ pc^ 

where is the observing frequency in GHz. Solving for SAV in (5) and (6), and equating, 
yields: 

(7) Leo = 3 X 10- VL'c,o Lq 
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The cubic dependence on Vr comes from the fact that L'^qi being in K, is independent of 
transition (for constant brightness temperature, T^), hence the luminosity wih increase as 
due to the definition of Tg, and another factor of Vr due to the increased hne width in Hz 
for the higher transitions, for a fixed velocity width in km s~^. 

Substituting for L'qq in equation (7) using (4) then yields: 

(8) Lco^^^l^'^^'ylLFiR Lo 

2.3. The CO emission from galaxies that reionize the Universe 

Again, making the assumption that the individual galaxy relationships apply to the 
volumetric average of galaxies on large scales, we can combine equations 8 and 3 to obtain 
the CO luminosity per unit volume required to reionize the neutral IGM: 

(9) pz.co=6.4xlO-V/e.c(l + ^)'(<^/5) LqMpc-^ 

Equation (9) is the fundamental relationship for the expected CO luminosity per co- 
moving cosmic volume dictated by the star formation rate density required to reionize the 
neutral IGM, based on empirical relationships between Leo, Lpm, and star formation rate 
for star forming galaxies at low and high redshift. 

We can go further, and consider observable quantities. Equation (5) gives the rela- 
tionship between observables S and Ay and the CO luminosity for a given galaxy. Again, 
assuming this applies to a volumetric average of galaxies (ie. Plco)i using equation 
(9) and (5), we obtain a velocity integrated CO flux per unit comoving volume required for 
reionization: 

(10) psAv = 6.4 X lQ-^\lD-/f-l{C/b) Jy km s"^ Mpc"^ 
where Da is the angular diameter distance in Gpc. Da ~ 1-0 ± 0.2 Gpc for 2; = 6 to 10. 

3. Examples 

Let us use equation (10) in an example. Consider a survey at 5' resolution at ^ = 8 for 
the CO 2-1 line (230GHz rest frequency, 25.6GHz observing frequency). The angular size 
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corresponds to 13 Mpc (comoving), which imphes a = 0.044 for the standard cosmological 
expansion. The comoving cosmic volume covered is then 2200 Mpc^. Multiplying through 
by this volume, equation (10) then yields SAV = 0.073(/esc/0.1)^^(C/5) Jy km s~^. The 
Az due to cosmic expansion corresponds to a depth of 1500 km s~^, yielding a mean signal 
over the band of: S = 50{fesc/0-l)~\C/5) /xJy. 

We can then use the standard Rayleigh- Jeans relationship to derive observed brightness 
temperature: Tg = 13605* X^6~'^ K, where 6 is the angular size in arcseconds, S is the flux 
density in Jy, and is the observing wavelength in cm. For a 50/iJy signal at 5' resolution 
and observing at 1.2cm wavelength: Tg = l.l(/esc/0.1)^^(C/5) /xK. 

Performing the same calculation at 15' resolution and Az = 0.13, yields: S = 450(/esc/0.1)~"^(C/5) 
/iJy, and again, Tg = l.l(/esc/0.1)^^(C/5) /iK. Likewise, for CO 1-0 at 15' resolution, the 
values are S = 110(/esc/0.1)-^(C/5) /xJy, and again, l.l(/e,c/0.1)-HC/5) /^K. 

Based on this calculation, we can estimate the expected mean CO brightness tempera- 
ture (observed-frame) as a function of redshift, for galaxies that can reionize the neutral IGM 
at a given redshift. The result is shown in Figure 1. Note that this is independent of tran- 
sition (assuming constant Tb for the low order transitions), and independent of resolution 
(being a mean surface brightness). 

4. Discussion 

I have considered the requisite CO emission coming from the galaxies that reionize the 
Universe at a given redshift. The calculation above leads to a mean surface brightness at 
z = 8oi: 

with a gradual variation with redshift (Figure 1). 

Note that real observations will entail detecting fluctuations in the CO signal around 
the mean. Also, this value of Tg is dictated by the star formation rate density at a given 
redshift that is needed to reionize, and keep ionized, the neutral IGM (equation (10)). In 
reality, the star formation rate density will evolve as it will, and we will observer over what 
redshift range the IGM reionizes. 

There are clearly substantial uncertainties in the calculation above. Explicit uncertain- 
ties involve /esc and C, both of which have current estimates that could range by a factor of 
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a few (see section 2.1). Determining the escape fraction of ionizing photons will be a major 
goal of future observations of z > 7 galaxies with JWST, while determining the clumping 
factor is clearly a goal of next generation HI 21cm reionization experiments. 

Implicit to the calculation are the relationships between star formation rate, FIR lu- 
minosity, and CO luminosity. The dust and CO production and heating in the first galax- 
ies remains uncertain. Hence, it remains unclear whether one can extrapolate standard 
Kennicutt-Schmidt star formation laws, and FIR luminosity to star formation rate relations, 
to very early galaxies. Detection of strong CO line emission, and thermal emission from 
warm dust, from z > 6 quasar host galaxies suggests that CO and dust can be produced 
rapidly in the first star forming galaxies, although these systems are at the high mass end of 
the galaxy distribution (Wang et al. 2008, 2010). As for CO excitation, clearly the CMB will 
'depopulate' the lowest levels at high z, although a calculation of the Boltzmann distribution 
suggests only a factor two or so decrease in the population of the first excited state at z ~ 7 
in a typical star forming galaxy, relative to having no background at all. Calibrating the 
relationships between CO luminosity, FIR emission, and star formation from z > 7 galax- 
ies will be a major goal with ALMA, the EVLA, and the JWST, in the coming years, on 
a galaxy-by-galaxy basis. This information can then be incorporated into the calculation 
above for a better estimate of the average CO emission seen on very large scales by intensity 
mapping. 

Another difficulty will be line confusion. The observed line cubes will contain different 
transitions at different redshifts. To make full use of the data a dual- frequency experiment 
may be required, probing two transitions of CO to eliminate ambiguities, although cross 
correlation with the HI 21cm images will also help. 

Overall, the estimate of Tg for CO above should be considered at best accurate to an 
order of magnitude. We can compare this value to more involved calculations in the recent 
literature. Righi et al. (2008) have considered the background radiation expected from CO 
lines from galaxies from z — 1 to 10, in the frequency range of ~ 10 to 50 GHz. Their 
calculation adopted the extended Press-Schecter formahsm to set the galaxy merger rate, 
and used analytic recipes for the implied CO emission from these galaxies. They obtain a 
mean signal from a given CO transition at z ~ 2 within a factor two of the mean calculated 
herein. More recently. Gong et al. (2011) have used the Millennium simulation results of 
Obreschkow et al (2009) to estimate the CO surface brightness fluctuations on large scales 
at z > 6. The Obreschkow et al. study involved generating mock galaxy catalogs in CO 
emission based on the Millennium dark matter cosmological simulation, and again, invoking 
recipes to extrapolate to CO luminosity for a given halo. Gong et al. find a mean CO 
surface brightness of ~ 0.1//K to 0.7 /iK at z ~ 7 (depending on model assumptions), with 
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fluctuations a factor of a few lower than this on scales ~ 10' (see also Purlanetto et al. 2011, 
in prep). The roughly factor two agreement between these sophisticated calculations and 
the simple estimate presented herein is encouraging, given the very different approaches. 

Is it plausible to detect such a signal? Without going into detailed instrument design, 
we can consider receiver noise and brightness temperature (ie. the radiometry equation). 
For a 20K receiver system in 1000 hours and a 130 MHz channel (= 1500 km s~^ at 26GHz): 
5T = 20K/[3.6 X 10*^ x 1.3 x 10^^^ = 0.9 /xK. Hence, achieving a la brightness sensitivity 
is plausible in a long integration. 

A key point to keep in mind is that the initial studies will be statistical in nature, eg. a 
cross correlation with the HI 21cm images, or study of the autocorrelation of the CO cubes 
themselves. The challenge will be obtaining the combination of very wide field (> 100 dcg^) 
with fiK surface brightness sensitivity at resolutions ~ 1' to 10', at frequencies between 15 
and 45 GHz. Gong et al. (2011) consider this point in more detail. They estimate that 
a 1000 element array of 0.7m diameter antennas with a 1 GHz bandwidth and 30 MHz 
channels operating at 15GHz could detect the CO power-spectrum to high significance in 
3000 hours on comoving scales of k ~ 0.1 Mpc~^ to 1 Mpc"^. Bowman et al. 2011 (in 
prep) are also considering the design for such an experiment, including the relative benefits 
of interferometers versus focal plane arrays on single dish telescopes. 

Ultimately, the critically unique analysis will be the cross correlation with the HI 21cm 
image cubes (Gong et al. 2011; Furlanetto et al. 2011), which is well beyond the scope of 
this letter. My main conclusion is that the simple calculation presented herein, based on 

the star formation rate required to reionize the Universe, supports the more sophisticated 
treatments of the problem, and suggests that the CO intensity mapping experiment at ^ > 6 
is plausible. 

CC thanks the Keck Institute for Space Studies for instigating this work, the CO inten- 
sity mapping design team for lively discussions, and the referee for insightful comments that 
improved this paper. 
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Fig. 1. — The expected mean CO brightness temperature (observed- frame) for galaxies that 
can reionize, and keep ionized, the neutral IGM at a given redshift, based on equation (10) 
and assuming fggc = 0.1 and C — 5. 



